Association of -31T>C and -511 C>T polymorphisms in the interleukin 1 beta (IL1B) promoter in Korean keratoconus patients.

Purpose To investigate the genetic association between unrelated Korean keratoconus patients and interleukin 1 alpha (IL1A), interleukin 1 beta (IL1B), and IL1 receptor antagonist (IL1RN) gene polymorphisms. Methods We investigated the association between IL1A (rs1800587, rs2071376, and rs17561), IL1B (rs1143627, rs16944, rs1143634, and rs1143633), and IL1RN (rs419598, rs423904, rs424078, and rs315952, variable number tandem repeat [VNTR]) polymorphisms in 100 unrelated Korean keratoconus patients. One hundred control individuals without any corneal disease were selected from the general population. Polymerase chain reaction (PCR) – restriction fragment length polymorphism (RFLP) analysis and direct sequencing were used to screen for genetic variations in the IL1 gene cluster. Haplotypes for the IL1 gene cluster were constructed using Haploview version 4.0. Results We analyzed a total of 12 polymorphic sites in the IL1 gene cluster. Among them, the −511 (rs16944) and −31 (rs1143627) positions in the promoter region of IL1B were significantly different between patient and control groups. The C allele of rs16944 (−511C>T, p=0.022, odds ratio of risk [OR]=1.46, 95% confidence intervals [CI] 0.94<2.27) and the T allele of rs1143627 (−31T>C, p=0.025, OR=1.43, 95% CI 0.92<2.22) were associated with a significantly increased risk of keratoconus in Korean patients. Linkage of the two alleles, −31*C and −511*T, was associated with an increased risk for keratoconus with OR=2.38 (p=0.012, 95% CI=1.116–5.046). The *C/*A genotype of rs2071376 in IL1A intron 6 was significantly different between the keratoconus patients and control subjects (p=0.034, OR=0.59, 95% CI 0.32<1.11). Other polymorphisms did not show an association with keratoconus risk. Conclusions This is the first report of IL1 gene cluster mutation screening in Korean keratoconus patients. Significant differences in allelic frequency of IL1B between keratoconus patients and the control group suggest that IL1B polymorphisms may play a role in the susceptibility of unrelated Koreans to develop keratoconus.

IL1 is a pleiotropic cytokine. It is involved in the inflammatory response, cell growth, and tissue repair in the cortex. The IL1 superfamily consists of three members, IL1 alpha (IL1α), IL1 beta (IL1β, predominant form) and IL1 receptor antagonist (IL1Ra), which are encoded by IL1A, IL1B, and IL1RN, respectively [17,18]. The IL1 superfamily genes are located in tandem in a cluster on chromosome 2q14. This gene cluster contains several polymorphisms [19]. Among them, some polymorphisms are located within the regulatory regions of the genes. Their localization in regulatory regions suggests that they may modulate IL1 protein production by directly affecting transcription, leading to their association with altered levels of IL1 [18,[20][21][22].
Biallelic polymorphisms in the IL1 gene cluster were determined by PCR-restriction fragment length polymorphism (RFLP) analysis and direct sequencing ( Table  1). The IL1RN VNTR polymorphism was evaluated by identifying the number of repeats in the (86)n sequence using agarose gel electrophoresis and direct sequencing [19,24,25].
The Hardy-Weinberg equilibrium (HWE) was calculated using the GenePop web version 4.0 program. To determine statistically significant differences in genotype and allele frequencies between the two groups, we used the χ 2 test or Fisher's exact test for the 2×2 contingency table file. The descriptive statistics for observed differences in allele or genotype distribution with the corresponding p values were analyzed using the JavaStat web software in combination with StatXact-8 software (Cytel Inc., Cambridge, MA). The strength of the association was estimated by odds ratio of risk (OR) and 95% confidence intervals (CI). Haplotype frequencies and associations were calculated with Haploview (version 4.0) that uses the expectation maximization (EM) algorithm [26]. Haplotype distributions were evaluated by the permutation test on the basis of 10,000 replications to obtain the empirical significance [26]. Values of p<0.05 were considered statistically significant.

RESULTS
We analyzed 11 SNPs and one VNTR marker in the IL1 gene cluster: rs1800587, rs2071376, and rs17561 for IL1A; rs1143627, rs16944, rs1143634, and rs1143633 for IL1B; and rs419598, rs423904, rs424078, rs315952, and the (86)n VNTR marker for IL1RN ( Figure 1). The IL1A, IL1B, and IL1RN genotypic and allelic frequencies in keratoconus patients are listed in Table 2 and Table 3. The genotype distributions of the variants of the 12 polymorphisms in the Assay method a : Screening method for polymorphic site and lists of restriction enzyme; Agarose: Agarose gel electrophoresis; Genotypes b : It was indicated restricted sizes by RFLP analysis and amplified products for VNTR.
IL1 gene cluster among the control subjects and the keratoconus patients were in Hardy-Weinberg equilibrium.
Of the four SNPs in IL1B, the rs1143627 (−31T>C) and rs16944 (−511C>T) positions in the promoter region of IL1B  Of the three SNPs in IL1A, the *C/*A genotype for rs2071376 in intron 6 was significantly different between keratoconus patients and control subjects (p=0.034, OR=0.59, 95% CI 0.32<1.11), but the allele frequencies were not significantly different between the groups. Other SNPs and the VNTR showed no significant genotype or allele frequency differences between keratoconus patients and control subjects.
In the haplotype analysis, we identified three haplotype block sets, rs2071376-rs17561-rs1800587-rs1143634 for Block A, rs1143627-rs16944 for Block B, and rs419598-rs423904-rs454078 for Block C (Figure 2). The T-C haplotype in Block B was more prevalent in keratoconus patients than in control subjects and was associated with a higher risk of developing keratoconus (p=0.034, OR=1.58, 95% CI 1.018<2.448). In contrast, the T-T haplotype in Block B was less frequent in the keratoconus group compared with the control group (p=0.009; Table 4).

DISCUSSION
Keratoconus is an abnormality of the interactive epithelialstromal system that causes a shift in the delicate balance between keratocyte proliferation and apoptosis [12]. Using the TUNEL assay, Kim et al. [12] confirmed the presence of keratocyte apoptosis in keratoconus corneas and the absence of keratocyte apoptosis in normal corneas. Keratocyte apoptosis may be triggered by increased basal IL1 release [12,14,27]. Previous studies indicated that corneal fibroblasts from keratoconus patients express about three times more IL1A mRNA than those from a normal cornea [15,28]. Zhou et al. [29] reported that the IL1 protein level was higher in the epithelium and endothelium of the keratoconus corneas than in normal corneas. In addition, IL1 upregulates keratocyte expression of collagenases, metalloproteinases, and other enzymes [30,31]. These enzymes have important roles in the remodeling of collagen during corneal wound healing. For example, persistent rubbing of the eye is likely to produce corneal epithelial trauma and increase the release of IL1 [15,16,32,33]. This would provide a unifying explanation for the association between keratoconus and factors associated with mechanical corneal epithelial cell injury such as poorly fitted rigid contact lenses, excessive eye rubbing, and allergic ocular surface disease [2,4,5,34]. Therefore, it has been suggested that interleukins are central modulators of the response to corneal injury.
The IL1 superfamily includes three secreted glycoproteins, IL1α, IL1β, and IL1Ra. Two of these, IL1α and IL1β, are biologically active while the third member, IL1Ra, is a receptor antagonist that modulates the effects of IL1α and IL1β [18,19]. IL1α and IL1β proteins are synthesized by a  variety of cell types including activated macrophages, stimulated B lymphocytes, fibroblasts, and endothelial cells.
They are powerful mediators of inflammation and the immune response. IL1A,IL1B, and IL1RN map to human chromosome 2q14. Gene variations that alter gene function are more likely to influence phenotypic characteristics such as risk of disease [18]. An in vivo study showed that the IL1A-889 *T allele was associated with increased IL1α and IL1β protein levels. Furthermore, ex vivo analysis of lipopolysaccharide (LPS)stimulated peripheral blood mononuclear cells indicated that production of IL1α from IL1A with the IL1A-889*T allele increased [35,36]. The two promoter variants of IL1B, −511C>T and −31T>C SNPs, have been repeatedly associated with multiple clinical conditions [37,38] such as cardiovascular disease [39] and gastric cancer [40] as well as with clinically observed differences in the levels of IL1β protein in vivo [39,41]. The *T allele at −31 in the TATA box of the IL1B promoter region is suspected to enhance gene expression and induction of IL1β [42]. The -31 *C allele was associated with increased IL1β levels and decreased IL1α levels [43]. It has been reported that the T and C variants at the IL1B −31 position regulate gene expression and differential binding of proteins. A specific haplotype, which is composed of the T allele at −511 and the C allele at −31, was significantly associated with a twofold to threefold increase in LPSinduced IL1β protein secretion [21]. There was also nearly a complete linkage disequilibrium between −511C>T and −31T>C in IL1B based on a cis interaction [40,42]. Some studies have suggested that the allelic polymorphism located within intron 2 of IL1RN play a role in differential modulating IL1 activity. The IL1RA*2 is associated with increased IL1Ra protein production in vitro [44]. The *C allele of +8006 C/T in exon 2 of IL1RN is associated with the VNTR allele 2 and is associated with lower expression of IL1Ra [45,46]. The IL1B-511C, IL1B-31T, and ILRN+8006T haplotype (C-T-T), which has been associated with lower levels of IL1β expression [21,38], was present in the extended protective haplotype. This extended protective haplotype contains the IL1RN+8006T allele, which is associated with elevated IL1Ra expression [45]. Therefore, these studies suggest that IL1 gene cluster polymorphisms may affect susceptibility to the development of specific diseases and that linkage disequilibrium of functional polymorphisms may indicate increased risk of developing these diseases.
In the present study, we have reported the results of a mutation screening of interleukin 1 genes in unrelated Korean patients with keratoconus. In the screen of four SNPs of IL1B, the genotypes of -511 *C/*C and -31 *T/*T of IL1B SNPs were associated with a 1.96 fold and 1.92 fold greater risk of developing keratoconus, respectively, (p=0.025 and p=0.027, respectively), showing a strong statistical association with keratoconus in unrelated patients. Two IL1B SNPs, −31C/T and −511T/C, have been implicated as potential risk factors for keratoconus among unrelated patients. Of the three SNPs in IL1A, the *C/*A genotype of rs2071376 in intron 6 differed significantly between the keratoconus patients and control subjects (p=0.034, OR=0.59, 95% CI 0.32<1.11), but the allele frequencies were not significantly different between the groups.
The T-C haplotype of IL1B −31 and -511 was more prevalent in keratoconus patients than in control subjects and carried a higher risk of keratoconus (p=0.034, OR=1.58, 95% CI 1.018<2.448). The T-T haplotype of IL1B −31 and −511 was less frequent in keratoconus patients than in the control individuals (p=0.009), which was associated with protection.
Our study is the first to report on IL1 gene cluster polymorphisms in Korean keratoconus patients. We detected two disease-associated SNPs, rs1143627 (−31T>C) and rs16944 (−511C>T), in IL1B and a disease-associated haplotype, the IL1B-31T/-511C haplotype, which is associated with regulation of IL1α and IL1β production. It proves that IL1 may have important roles as central modulators of keratoconus development. In conclusion, the present study suggests that two SNPs in IL1B predict keratoconus predisposition in unrelated Korean patients, although further research is necessary to elucidate the relationship between the expression levels of IL1α, IL1β, and IL1Ra and SNPs of IL1B in keratoconus patients.